Kink dynamics in a one-dimensional growing surface 
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A high-symmetry crystal surface may undergo a kinetic in- 
stability during the growth, such that its late stage evolution 
resembles a phase separation process. This parallel is rigorous 
in one dimension, if the conserved surface current is derivable 
from a free energy. We study the problem in presence of a 
physically relevant term breaking the up-down symmetry of 
the surface and which can not be derived from a free energy. 
Following the treatment introduced by Kawasaki and Ohta 
[Physica 116A, 573 (1982)] for the symmetric case, we are 
able to translate the problem of the surface evolution into a 
problem of nonlinear dynamics of kinks (domain walls). Be- 
cause of the break of symmetry, two different classes (A and 
B) of kinks appear and their analytical form is derived. The 
effect of the adding term is to shrink a kink A and to widen 
the neighbouring kink B, in such a way that the product 
of their widths keeps constant. Concerning the dynamics, 
this implies that kinks A move much faster than kinks B. 
Since the kink profiles approach exponentially the asymptot- 
ical values, the time dependence of the average distance L(t) 
between kinks does not change: L(t) ~ lnt in absence of 
noise, and L(t) ~ t 1 ' 3 in presence of (shot) noise. However, 
the cross-over time between the first and the second regime 
may increase even of some orders of magnitude. Finally, our 
results show that kinks A may be so narrow that their width 
is comparable to the lattice constant: in this case, they in- 
deed represent a discontinuity of the surface slope, that is an 
angular point, and a different approach to coarsening should 
be used. 



I. INTRODUCTION 

A crystalline surface growing under a flux of incoming 
particles from the vapour phase represents a typical ex- 
ample of an out-of-equilibrium system. Its microscopic 
evolution may be described as follows: once the adatom 
has arrived on the surface, it performs a thermally ac- 
tivated diffusion process till it is "trapped" somewhere, 
or it evaporates by coming back to the vapour phase. 
Which surface relaxation mechanism (surface diffusion or 
evaporation/condensation) indeed prevails, it depends on 
the temperature and the specific parameters of the ma- 
terial jjj . Anyway, for a wide class of materials (mainly 
metals), at the relevant temperatures for Molecular Beam 
Epitaxy (MBE) desorption may be neglected. In this 
case, two different "traps" may be effective: another 
adatom (giving rise to a nucleation phenomenon), or a 
step. If we limit ourselves to the case of a high-symmetry 
surface, there are no preexisting steps and therefore the 
mentioned step belongs to a growing island. Once islands 



have coalesced, leading to the completition of one layer, 
steps should disappear and the whole previous process 
should start again. 

The previous qualitative picture applies to the case of 
a stable layer-by-layer growth. In reality, it is hindered 
both by noise and by possible instabilities; sources of 
noise are fluctuations in the flux of incoming particles 
(shot noise), in the surface diffusion current (diffusion 
noise), and in the nucleation events (nucleation noise). 
Whilst the first two have been well-studied in the con- 
text of several different models B, the latter one still 
needs a more basic comprehension J3| . Concerning deter- 
ministic instabilities, the main responsible —and perhaps 
the sole, for a homoepitaxial high-symmetry surface— of 
the destabilization of the flat surface, is now known as 
Ehrlich-Schwoebel effect [Q: an adatom approaching a 
step from above or below may have different probabili- 
ties of attachement. If the sticking from above is discour- 
aged, an adatom which has the possibility to choose be- 
tween two different kinds of steps (an ascending one and 
a descending one) will stick preferably to the ascending 
one, thus determining an up-hill current. It is important 
to remark that this is a purely out-of-equilibrium effect, 
because at equilibrium detailed balance forbides such a 
current! 

Even without entering into details, as it will be done 
in the next section, it is possible to explain here the ef- 
fect of such a mechanism. In fact, as was firstly pointed 
out by Villain j^j, the resulting surface current j = vm 
(m = dz/dx being the local slope of the surface and 
z the local height), once put in the evolution equation: 
dz/dt = —dj/dx gives rise to a diffusion- type equation: 
dtZ = —vd%z, where the negative sign of the diffusion 
constant {—v) is responsible for the instability of the flat 
surface (z = constant). In the present paper, we will 
mainly be concerned with the late stages of this insta- 
bility, when additional and nonlinear terms must be in- 
troduced to describe the dynamics of the surface. In the 
next section, we will introduce a more general expression 
for the surface current j and we will take into account 
the breaking of the z — > — z symmetry, induced by the 
flux F of atoms. 



II. THE SURFACE CURRENT 

The study of a growth process may ideally be divided 
into two main steps: the first one starts from some mi- 
croscopic point of view and should arrive to a contin- 
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uum description of the surface/interface; the second one 
may even assume a given evolution equation and study 
it. Most of the difficulties encountered in a theoretical 
study of MBE are related to the first step. In the present 
paper, we will limit ourselves to a one-dimensional high- 
symmetry surface, and in this section we will introduce 
and justify a specific Langevin-type equation. 

The local surface height z(x,t) is generally supposed 
to satisfy a local equation of the form dtz = aF + 
!F{d x z, O^Zi ■ ■ ■), F being the incoming flux, a the in- 
plane lattice constant, and T a function of the local pro- 
file of the surface (the out-of-plane lattice constant is put 
equal to one, i.e. z is adimcnsional) . The underlying hy- 
potheses have been discussed in Ref. Q, where we have 
shown that the appearance of angular points in the sur- 
face profile may be treated correctly solely through the 
introduction of a nonlocal equation. We will take up this 
point again, at the end of the article. 

The flux F contains a constant part Fq, which is "elim- 
inated" by redefining z(x,t): z —* z — aF$t, and a fluc- 
tuating part 5F(x, t) which represents the so-called shot 
noise, which is supposed to follow a Gaussian-like distri- 
bution: 

(SF(x,t))=0 (1) 
(SF(x, t)SF(x', t')) = 2F S{x - x')5{t - t') (2) 

In the limit of negligible desorption, and if overhangs 
are forbidden, surface growth proceeds by conserving 
both mass and volume: therefore, the function T must 
be derivable from a surface current j, and the evolution 
equation will be written in the form: 



dtz(x, t) = —ad x j + a5F(x, t) 



(3) 



The central question is which current j governs the 
evolution of the surface: a still debated question even 
for the simplified model of a one dimensional surface, 
as shown by the following discussion on the different 
terms appearing in j. Symmetry arguments simply tell 
that j does not depend on z |j, but on its derivatives 
(to = d x z, m! = d\z, to" — d^z, . . .) and that —on a high- 
symmetry surface— it must be an odd function of x: so, 
a term proportional to m or to" satisfies this request, but 
if proportional to m', it does not. 



A. Ehrlich-Schwoebel current 

In the Introduction, we mentioned the Ehrlich- 
Schwoebel effect, which gives rise to a slope-dependent 
current: j_es(w). Since it must be an odd function of to, 
its form at small to will be jes = vm - The coefficient v 
depends HH on the flux F , the diffusion length to and 
the Schwoebel length io measures the typical linear 
distance travelled by the adatom before meeting another 
one and forming the nucleus of a growing island. It rep- 
resents the "maximal" size of a terrace, because if i > Id 



the probability to nucleate a new island on it is very high; 
during the first stages of growth, when the surface is still 
more or less flat, Id is also the typical size of a terrace. 
It is not so when an instability develops: in this case, I 
may be much smaller than Id and the slope m = 1/1 
may be fairly large. Indeed, the slope 1/£d discriminates 
between a nucleation-dominated regime (m <C and 
a step- flow regime (m 1/£d)- the latter is generally 
relevant for vicinal surfaces, which grow through sticking 
of adatoms to preexisting steps; anyway, if a flat surface 
develops an instability with regions of high slope, such 
regime becomes important also for high-symmetry orien- 
tations. 

The second relevant length, the Schwoebel length, is 
a measure of the asymmetry in the sticking coefficients 
of an adatom to a step. Its simplest form |jj is: is = 
a(D-/D+ — 1), where D + and D_ are such coefficients 
for an adatom approaching the step from above (-D+) 
and below (D-). The existence of an Ehrlich-Schwoebel 
effect means that D + < £>_ and therefore is > 0. To 
describe the meaning of is, let us consider a terrace of 
size I (< £d): if is *C i, only a fraction is/i of the 
fallen adatoms will contribute to the uphill current, and 
therefore jes = {F i)(is/i) = F is, since the number 
of atoms arriving per unit time on a terrace of size i is 
nothing but Fgi. Conversely, if is 3> i a U the adatoms 
will stick to the ascending step, and so: jes — FqI. The 
simplest interpolation formula, valid for any value of is 
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F is 



l + is\m\ 



H > l/h 



(4) 



This formula also allows to obtain a semiquantitative ex- 
pression for the parameter v. in fact, when \m\ ~ 
Eq. (Q) must match the expression valid at small slopes: 
Jes = vm. The result is: v = F i s i D /(is + #d)- It is 
important to remark that all the previous considerations 
may be made more rigorous but in this section wc 
are mainly interested in justifying the expression for the 
current j, rather than in deriving it. 

The main characteristic of the Ehrlich-Schwoebel cur- 
rent just discussed is that it has no zeros other than 
m = and m = ±oo. A zero in jes is extremely im- 
portant |[o|], because the other terms in j will be seen 
to depend on higher order derivatives of z(x,t). So, a 
constant slope mo may be a stationary slope if and only 
if J-Es(mo) = 0. An extra-zero tuq may have different 
origins: the symmetry of the crystal lattice |lo|Jll|| , non- 
thermal relaxation mechanisms or a transient mo- 
bility of the adatom just after the deposition For 
example, the slope at 45 degrees corresponds in a cubic 
lattice to the high- symmetry orientation (11): we expect 
that jes vanishes on it, as it vanishes on the (10) (to = 0) 
and (01) (to = oo) orientations. A different example is 
the following: if atoms falling in the vicinity of a step 
have a higher probability to land on the lower terrace, 
or to kick down the step adatom, a down-hill current j± , 
proportional to the density of steps and therefore to the 
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slope to, will appear: j± = —v'm. So, if v > z/, a zero 
will appear when (Jes + 3l) — 0. 

Whatever is the origin of extra-zero(s) in the slope- 
dependent current, we can introduce two different mod- 
els, according to the presence (model I) or absence 
(model II) of zeros at finite slope. The simplest expres- 
sions of jes for the two models, having the correct sym- 
metry properties are [13]: 



model I 
model II 



jes = vm(l - to /m ) 
vm 

Jes 



1 



(5) 
(6) 



D 



Model II does not correspond to a phase sepa ration pro- 
cess (see Sec. Ill): it will be discussed in Sec. VII]. 



B. Mullins-like current 

The most famous "equilibrium" current is perhaps the 
one (j'm) introduced by Mullins |HJ forty years ago, 
to study the relaxation towards equilibrium of a non- 
singular grooved surface. A simple derivation starts from 
writing jm as the gradient of a chemical (surface) poten- 
tial: jm = — r^jcju, where Y is the adatom mobility, and 
afterwards to derive fi from a surface free-energy: 



S£ 

5z(x) ' 



with £ = a J dx\J I + m 2 a 2 (7) 



By combining the different equations, in the limit of small 
slopes we obtain: 



j M = Km"(x) 



(8) 



with K = a 2 Ycr. 

The usage of this expression in our problem may be 
questionable in at least two respects: first, it applies to 
a nonsingular surface, i.e. above the roughening transi- 
tion Tr; second, it applies to a close-to-equilibrium sur- 
face. Concerning the first remark, our surface is a high- 
symmetry one and therefore almost necessarily below Tr, 
because for a high symmetry orientation the roughen- 
ing temperature is equal or nearly equal to the melting 
temperature T^, while ordinary temperatures for MBE 
are well below Tm- Nevertheless, our surface —which is 
strongly out-of-equilibrium— contains a lot of steps be- 
cause the incoming flux makes the surface rough [ jl6| : 
therefore, the surface current should be nonsingular at 
zero slope. 

The latter remark is more "critical" : the Mullins cur- 
rent derives from thermal detachment of atoms from 
steps, in order to minimize the surface free-energy. It is 
not clear if such process is effective in presence of a flux 
F. For example, Stroscio and Pierce |17j state that ther- 
mal detachment is negligible in the homoepitaxial growth 
of Fe (at least at room temperature) and therefore they 
do not write |l§| ] such a term in the current. Anyway, 
it has been shown Jqjlq] that the current (||) may derive 



also from nonequilibrium effects: nucleation noise and 
diffusion noise. The first one should be dominant and 



correspond to the value |l 16 1: K — F l 4 D . 



C. Symmetry-breaking current 

The terms in the surface current which have been intro- 
duced so far not only satisfy the x — > — x symmetry (be- 
cause j(—x) = —j(x)), but they also fulfil the up-down 
symmetry, corresponding to the change of sign of z. In 
fact, if z — > —z both jes and jm change sign. However, 
there is no reason to expect that surface growth proceeds 
by conserving such symmetry, since the flux breaks it. 

A symmetry-breaking (SB) term is intrinsically non- 
linear, because any current of the form j ~ d"z(x,t) 
changes sign with z. The lowest order expression which 
changes sign with x but does not change sign with z is: 



jsB = d x A(m 2 ) 



(9) 



where A is any even function of the local slope. The 
simplest form for A: A — (A/2)m 2 has been introduced 
by Sun et al. JlSj] . It is also called "conserved Kardar- 
Parisi-Zhang term" , because in Eq. (||) it looks like the 
laplacian of (d x z) 2 , i.e. the nonlinear term of the KPZ 



equation |20 |. 

The current (|9|) is not derivable from a free energy. 
As pointed out by Somfai and Sander |pl[ it is necessary 
to rise the order of jsB to make it derivable from some 
free energy (for example, j SB ~ d x [(m') 2 } = j dxT sb 

with T SB ~ K) 3 )- 

Before proceeding, let us discuss the physical origin 
°f 3 SB ■ When there is a gradient in the density p of 
adatoms, a current of the form j = —Dd x p is expected, 
where D is the diffusion constant. In the case of a growing 
surface, the applicability of the previous expression is not 
obvious, because steps are sink for diffusing atoms and 
—at least if thermal detachment is forbidden— interlayer 
diffusion is absent. In spite of this, the above expression 
may help in understanding: in fact, adatom density on 
a terrace depends on its size £, because a larger terrace 
collects more atoms from the flux than a smaller one. 
So, p — p{€) — p(\m\). In other words, the function A 
appearing in j$B seems to be proportional to the adatom 
density itself. 

This interpretation can be made more rigorous for 
large slopes (|to| = l/l > 1/Id), where nucleation of 
new terraces is absent and p can be simply determined 
by solving the diffusion equation dtp = Fq + Dd 2 p in the 
quasi-static approximation [dtp — 0) and with p(0) = 
p(£) = as boundary conditions (i.e. steps are perfect 
sinks). The resulting average density on the terrace is p ~ 
{F a /D)( 2 and the current is j SB ^ -F Q d x (l/m 2 ). This 
expression agrees with those determined, with different 
methods, by Politi and Villain || and by Krug |2^] . Hunt 
et al. suggest that jsB may derive from the sticking- 
asymmetry induced by the Ehrlich-Schwoebcl effect: nev- 
ertheless, jsB does not vanish even if Is — [pLp2|- 
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One could ask why the average value of p is taken. The 
answer is that inhomogeneities in the adatom density on 
a given terrace give rise to jes* I n fact, if no Ehrlich- 
Schwoebel effect is present, p(x) is symmetric with re- 
spect to the center of the terrace and therefore the av- 
erage value of d x p vanishes. Conversely, if is > then 
(d x p) terrace ^ anc ^ ^ corresponds just to jes- This 
remark stresses the "similar" origin of jes an d jsB- It 
is likely that a systematic derivation of the surface cur- 
rent should give all the terms we have introduced: ]es 
(which depends on the slope to), j$B (which depends on 
the curvature to'), and j'm (which depends on a higher 
order derivative: to"). Anyway, a rigorous derivation is 
still lacking at the moment, above all for a high-symmetry 
orientation. 



D. The current of our model 

In the following, we will study the dynamical evolution 
of the surface, as determined by the current: 



3 = 3es + 3 m + 3 SB 



where: 



Jes 



1- — 



j M = Km" 

]SB = Atoto' 



(10) 



(11) 

(12) 
(13) 



The reason of our choice is clear: we want to study 
the effect of the symmetry-breaking current (jg b ) on the 
phase separation process determined by the other two 
terms of the surface current (Jes + 3 m), an d at this aim 
we choose the simplest expression for jes —which must 
have a zero at a finite slope too— and for jsB —for which 
we take A[m 2 ) = Ato 2 /2. In the last section, we will 
discuss how the conclusions depend or not depend on the 
present choice. 



III. EVOLUTION IN ABSENCE OF THE 
SYMMETRY-BREAKING CURRENT 

In the "language" of surface growth, the evolution of 
the surface proceeds as follows: after a time t* an in- 
stability of the flat surface with a well-determined wave- 
length L* develops. In this linear regime, L* is constant 
and the amplitude increases exponentially. Afterwards, 
because of the nonlinearity of Jes a coarsening process 
takes place: the wavelength L(t) of the mound- like (or 
pyramid-like) surface profile increases in time, whilst the 
maximal slope tends to the constant values ±too- So, the 
surface is "made up" of neighbouring regions where the 
slope is alternately (nearly) equal to +TOo and —too. 

The first stages of growth can be analyzed by lineariz- 
ing Eq. (||) with the current (|To|): 



a^dtzix, t) = ~vd 2 x z(x, t) - Kd A x z{x, t) (14) 

which shows |||)],^3| that the flat surface is unstable 
against deformations of wavelength larger than L c — 
2tt-v/ KJ v. The most unstable mode correspondes to 
L u = \/2L c and its amplitude grows as exp [(a^ 2 /4i^)t] . 
So, L* = L u and t* = (\Kjav 2 ). 

The nonlinear profiles of the mounds are determined as 
stationary solutions of Eq. (|J) , that is to say as solutions 
of the equation j = 0: 



3Es(m) + Km"{x) = 



(15) 



This equation can be derived by the following La- 
grangian: 

C = (K/2)m' 2 - V(m) with V'(m) = j E s{m) 

which corresponds to an anharmonic pendulum, once 
we have identified the slope to as its spatial coordi- 
nate and x as the time. Since the potential V(m) — 
(y /2)m 2 [l — to 2 /2toq] has two symmetric maxima in 
±too, the period of the oscillation (i.e. the wavelength of 
the surface profile) diverges when its amplitude (i.e. the 
maximal slope of the surface profile) goes to too. If Jes 
followed model II, V(m) would have no maxima and no 
limitation on the slope would be present. 

By going on with this mechanical analogy, the exis- 
tence of coarsening requires a condition on the station- 
ary configurations: the period of the oscillation must be 
an increasing function of the amplitude p4j ; a condition 
which is surely fulfilled by the potential V(ro), since the 
quartic correction has a negative sign! Clearly, coarsen- 
ing also requires that these stationary solutions are not 
stable: more precisely, they must be unstable with re- 
spect to wavelength fluctuations, but stable with respect 
to amplitude fluctuations. 

The previous mechanical analogy helps in understand- 
ing why the surface keeps a regular profile and also allows 
to determine this profile at a given time, but it is not ef- 
fective in determining the time dependence of L(t), i.e. 
the coarsening law Q. To this end, we must observe 
that the evolution equation for the local slope to (which 
represents the "order parameter" of our problem) satis- 
fies the noisy Cahn-Hilliard equation |25j|: 



' l dtm — d 2 



5T_ 

8m 



■r){x,t) 



where T 



dxC 



This equation corresponds to a phase separation 
process, where the order parameter is conserved 
(dt fdxm(x,t) = 0). The system is made up of do- 
mains where to equals one of the two degenerate min- 
ima of the potential energy [/(to) = —V(m); domains 
which are separated by domain walls move in order to 
minimize the "action" T . Domain wall (or "kink" ) move- 
ment is determined both by their (deterministic) interac- 
tion and by fluctuations induced by the conserved noise. 
We will see that the growing surface (even in presence of 
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the symmetry-breaking current jsb) can be mapped in 
a one-dimensional system of interacting kinks which an- 
nihilate, so that the average distance L(t) between kinks 
increases in time. 

By using this method for the symmetric case {jsB = 
0), Kawasaki and Ohta |2(| have found the equation of 
motion for the kinks, which has been then studied by 
Kawakatsu and Munakata [^tJ. The final result is that 
L(t) grows logarithmically with time if noise is absent 
and grows as i 1 / 3 if noise is present. 

IV. KINK PROFILES 

A stationary kink M(x) is defined as a monotonic so- 
lution of j[M{x)] = 0, with M(x) tending to (different) 
minima of U(m), when x — > ±00. In the present case, 
there are only two symmetric minima in ±mo and there- 
fore only two kinks M±(x) are possible, the subscript 
corresponding to the sign of its first derivative, i.e. to 
the curvature of the surface profile. 

The surprising result is that the "shape" of the kink 
does not change because of the introduction of the 
symmetry-breaking term. To see it, let us replace the 
expression: 

M±(x) = ±m t&nh(n±x/2) (16) 

in the differential equation j = 0: 

Km"(x) + um{l — m /m§) + Xmm = (17) 

We obtain the following second degree equation for the 
parameters n±: 

Kk± T \m K± - 2v = (18) 

which gives the positive solutions: 



K± = \^JX 2 m 2 + %vK ± Xm j /2K (19) 

Two limiting cases, corresponding to weak and strong 
symmetry-breaking, will be frequently used: 

Xm < \[%vK k + = k_ = \j2vjK = n (20) 
Xnio 3> v&vK k + = Xm^/K k_ = 2v/Xmo (21) 

So, the effect of jsB is two create two classes of kinks: 
kinks "A" , given by the profile M + (x) and characterized 
by a width (1/k + ), and kinks "B" , given by the profile 
M-(x) and whose width is (1/k_). For a strong jsB, 
k+ kinks A are much narrower than kinks B. It 

must also be observed that the product (k + k_) does not 
depend on A, since it equals (see the algebraic equation) 
{2v/K). In other terms, the effect of jsB is to shrink 
kinks A and to widen kinks B, in such a way that the 
product of their widths keeps constant. 



V. FROM SURFACE DYNAMICS TO KINK 
DYNAMICS 

In this section we will describe the method to solve the 
growth equation for the surface-slope profile: 

ar x d t m = D x [Km" - U'{m) + Xmm] B x = -d 2 x 

(22) 

in a "multi-kink" approximation. Since our approach 
follows that one introduced by Kawasaki and Ohta |26| ] 
to study the above equation in absence of the A-term, 
we will expose the main calculations in App. ^ and here 
we will limit ourselves to explain the general lines of the 
method. 

Once a kink is inserted in our problem, it moves with 
a given (constant) velocity v° and a profile m(x, t) — 
M(x — v°t), where v° is found by solving the eigenvalues 
problem obtained by putting m(x, t) in (|22|). Our system 
is made up of an ensemble of kinks A which alternate to 
kinks B, and we will look for an approximate solution of 
) as a superposition of kinks centered in Xi and mov- 
ing with velocity Vi. Because of the interaction between 
kinks, Vi is not a constant, and depends on the position 
of the other kinks. In principle, the nonlinear part of 
U(m) (i.e. the quartic term m 4 ) gives rise to terms of 
n-kinks interaction: we will adopt a "binary-interaction" 
approximation, which will be further simplified by limit- 
ing to nearest-neighbour interaction. This procedure is 
justified by the fact that we are interested in the late 
stages of growth, when the distance between kinks is 
much larger than the width of their cores (= 1/k±): so, 
they interact only through the tails of the profiles, which 
means that the interaction decays exponentially, since 
tanh(/c.x/2) ~ ±1 =F exp(— k|x|) when x — > ±00. For 
the same reason, the velocities Vi and the accelerations 
Vi will be considered "small" , because the typical size of 
the mounds grows slower than linearly: This means that 
the velocity of the coarsening process goes to zero, as 
time increases. 

As a final result, we obtain a Langevin equation for 
the discrete variables Xi(t), or — equivalently— for the 
kink- kink distances Xi(t) = Xi + \(t) — Xi(t), which will be 
studied by translating it in a Fokker-Planck equation. 

The treatment of Eq. (||) (see App. [§) gives the fol- 
lowing coupled equations for the kink positions: 

^a^ml^i-iy-^x, -Xj\xj = 
3 

(Ci) + (C 2 ) + (C 4 ) + Vi(t) (23) 

where: 

(d) = 8ism 2 [R p (X t ) - RpiXi-i.)] 

(C 2 ) = (3^/2,)ml^\[R^{Xi) - R^X^)] (24) 

(Ci) = -(3AmlK P X[R (Xi) - ^(AVi)] 

and: 
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(»*(*)) 







-AmlFoi-lf-^Xi-XjlSit-t') 



Let us explain the notations: The i-th kink is cen- 
tered in Xi, and —because of the breaking of symmetry- 
two different classes of kinks exist. In accordance 
with Sec. IV, their profiles are given by: Mp(x) = 
Pmo tanh(K^a;/2), where (3 = ±1. We will assume that 
the i-th kink is of class (3 (whatever is its value) and its 
nearest-neighbours of class —(3. The quantity: 



Rp(x) = exp(—K/3x) 



(25) 



in the (Cj) expresses the interaction between kinks, when 
the distances \xi±i — Xi\ are large compared to (l/np). 

Eq. ( |23| ) can also be written in matrix form: Aijij — 
Ii + r/i. The matrix A takes into account the kinematical 
coupling between kinks, due to the conservation of the 
order parameter, and / contains the forces between kinks. 
The matrix A can be inverted p7y , giving a tridiagonal 
and symmetric A^ 1 : 



A- 1 - 



1 



4ml V Xi 



4ml x i 



1 



The evaluation of A I is trivial: 



4m\ 



h + li+l . h + h- 



X, 



X, 



(26) 
(27) 

(28) 



and the explicit expression of Ii is found directly from 
Eqs. <M: 



Ii — Rg{Xi) — R%(Xi-i) 



(29) 



where R^(X) is a linear combination of the two different 

R/3(X): 

R*p(X) = c R (X) + dpR-p(X) with 



cp = &vml — f34m^np\ 



and 



d p = p(4/3)m 3 Kp\ 



Concerning the noise, it is preferable to work with 
quantities which are not spatially correlated. To this 
end, the matrix A -1 is written as the product PP T and 
new noise variables fj = P T r\ are defined. Since P is a 
bidiagonal matrix whose nonvanishing elements are: 



(30) 



rji is given by fji 
suits that: 



2mo \JXi 
\fa{r)i + r) i+ i)/(2m.oVXi), and it re 



(>*(*)) - (rjiitffjjit')) = 2aF 8 ij S(t - t') (31) 

In order to eliminate the constant factor in the correla- 
tor, we simply put: % = \/2aFo^i. This way, the final 
equation for kink dynamics is: 



±i(t) = 



4ml 



Ii + h+l . Ii + Ii-1 



X; 



2mQ 



Xi-i 



'Xi y/ Xi-1 



(32) 



VI. FORCES AND KINK VELOCITIES 

In this section we want to discuss the effect of the sym- 
metry breaking on the equations of motion for the kinks. 
Since we are here interested in the deterministic part of 
the interaction, we will not consider the noise. There- 



fore Ec 



takes the form: As 



L. Kawasaki and 



Ohta [E6[ suggest to look on Ii as the force acting on the 
i-th kink. Let us consider the two opposite limits: A = 
and Xm > VBvK. For R* p {X) we obtain: 

R%(x) = Svml exp(— kqx) 

in the first limit (A = 0), and: 

R* + (x) = {4\ 2 m^/3K)cxp{-n x x) 
R*L(x) = 16vmlexp(—K\x) 



in the second one (Atoo 3> V8vK) . In the previous equa- 
yj2v/K and k\ = (Xm /K) > k q . They 



tlOnS, Kq 

correspond to k_ in the two pertinent limits 

In the case of absence of the A-term, we simply get 



Ii = 8isml{exp(- n Xi) - exp(— kq^Q-i)] 



(33) 



This equation can be interpreted by saying that there is 
an attraction between kinks, proportional to exp(— kqX). 
If A 7^ (and "strong"), than we must distinguish be- 
tween positive and negative kinks: 



Ii = (4A 2 TO^/3K)[exp(- 
li = 16vml{exp(— K\Xi 



K X Xi) - exp(- 
— exp(— k\X, 



-nxXi-x)} 03 > 0) 
-i)] 09 <0) 



The first comment is that symmetry breaking implies 
that a positive kink is attracted (by a negative one) more 
strongly than a negative kink is attracted by a positive 
one. In other words, if we assign a mass to a kink, a 
negative kink weights more than a positive one, and the 
mass is proportional to the width of the kink itself. 

This interpretation seems to be satisfactory, but if we 
analyse the velocities Xi (t) rather than the "forces" Ii the 
picture becomes more complicated. In the limit A = 
we have: 



Xiit) = 2av 



exp(-K X t ) _ exp(-K Xi_i) 

Xi—i Xi 
exp(-K X i+ i) exp(-K Xi- 2 ) 



X; 



X 



i— 1 



(34) 



G 



So, the effect of the conservation law (i.d. of the ma- 
trix A) is that ii(t) depends not only on the positions 
of the nearest neighbour (nn) kinks (xi±i), but also 
on those of the next-nearest (nnn) ones (xi± 2 ). Even 
more important, the nnn- "interaction" is of the same or- 
der of magnitude as the nn-one! Whilst the interpreta- 
tion of exp(— KoXi) and exp(— KoXi-i) in Eq. (|33|), as 
—respectively— the interaction with the kinks + and 
(i — 1) is straightforward, in Eq. ( |34| ) the generic term 
exp(— kqXi) is divided by a different Xj, and therefore 
a similar interpretation becomes less evident. Anyway, 
if we don't ascribe too much importance to the quanti- 
ties Xj in the denominator, Eq. ( |34| ) says that kink i is 
attracted both by nn-kinks and nnn-kinks: the "interac- 
tion" between i and (i±2) has a kincmatical origin (con- 
servation of the order parameter) and indeed depends on 
(Xi±2 — Xi±i) rather than on {xi±2 — x{). A further com- 
ment is that in the evaluation of (ij + Zj+i) two terms 
cancel exactly, because in this case action and reaction 
are opposite and equal. 

If now we consider the case of a strong symmetry 
breaking term, the velocity takes the form: 



a\ 2 m,Q 
3K 



1 



1 ! 1 
Xi Xi-i 
1 



exp(-K A X l ) 



exp(-K A X i _i^ 



Xi Xi-i 

exp(-K A X l+ i) exp(-K A Ai_ 2 ) 



and: 

Xi(t)\fj <0 



-Aav 



a\ 2 rriQ 
3K 



X, 



Xi 



cxp(-K A Ai) exp(-K A X 4 _i) 



Xi 



Xt-i 



+ 



exp(-K A X i+ i) exp(-/i A Aj_ 2 ) 



Xi 



Xi-i 



(35) 



The surprising result is that the sign of the terms pro- 
portional to exp(— K\Xi) and exp(— K\Xi_i) is inverted: 
so —because of kinematics— a negative kink is subject 
to a repulsive interaction with its nn-kinks! This re- 
sult derives from the unbalancing of action and reac- 
tion. A closer inspection of the derivation of Eqs. ( |35| - 
|35|) allows to give the following interpretation: if /y 
means the force exerted by the kink j on the kink i (so 
that Ii — + /i,i-i), then kinematics determines 

that the effective force /j,i±i is a linear combination of 
(/m±i + fi±i,i) and fi,i±l- If A = 0, the first term van- 
ishes, but if A ^ it does not: furthermore, for a negative 
kink fi±x t i prevails over fi t i±\ and it corresponds to a re- 
pulsive force, for kink i. 

The conclusion we draw from the previous considera- 
tions is that negative kinks move much slower than pos- 
itive kinks. This results on one side from the fact that 
a bigger mass can be attributed to them, and on the 
other side that they are subject to an effective repulsive 
nn-interaction. 



VII. FROM KINK DYNAMICS TO 
COARSENING LAWS 

The interesting dynamical variables are the kink-kink 
distances X i; rather than the kink-positions X{. So, from 
Eq. (|32|) we obtain: 

**® = £l X 

' (B%(xi+ 2 ) + R*_ (x i+1 ) - R* fi {xi+i) - R*-p(xi)) 



Xi+i 
1 



Xi-i 



(R* (xi) + Rtpixi-i) - i$(xi_i) - Rtp(xi- 2 )) 



+- 



2mo 



6 



6- 



yX+i \/ Xi-i 



(36) 



The previous equations have the form: 

ft(t) - Ui({q}) + J2 &AUm with (37) 
3 

and we can therefore obtain a Fokker-Planck equation 
for the probability p({q}, t) of finding a given distribution 
{(?}, at time t. Two different procedures exist ^8|, due to 
Ito and to Stratonovich, but as remarked by Kawakatsu 
and Munakata [^7| the result is the same. This is true 
even in presence of the symmetry breaking term, because 
the two procedures may differ with respect of the term 
Gij, which does not change if the A-term is added. 
The Fokker-Planck equation writes: 

&i = " E TW^^ + 9 E E SkmStmP 



2 dq k dqe 



Its actual form, in our case, is 
dp \ - d 

AC 



[Uk{{X})p] 



(38) 



(39) 



Am 2 , 



E- 

4^ Xu 



d 2 



d 2 



o 2 



dx 2 _ x dx 2 +1 ax k ^dx k+1 



where Uk is nothing but the "deterministic" velocity of 
the fc-th kink. 

We are interested in the time dependence of the aver- 
age value of Xi (which does not depend on i). To this 
end, we define the distribution functions: 



/•OO 

g(X i ;t)= / (dX) lP 
Jo 

/■OO 

g 2 (X i ,X i+ v,t)= / (dXh^p 
Jo 

/>oo 

93(Xi,X i+1 ,X i+2 ;t) = / {dX) ll + ll + 2 p 
Jo 



(40) 
(41) 
(42) 
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The notation (dX) i l r j means that the integration is 
performed on all the variables Xj but Xi, X^i, .... 

The details of the calculation follow Ref. ]2^] and there- 
fore they will not be given here. By using the factoriza- 
tion approximation: 

g 2 (X h X i+1 ;t) = g(Xi)g(X i+1 ) (43) 
g 3 (X h X i+1 ,X i+2 ;t) = g(X t )g(X l+1 )g(X l+2 ) (44) 

and integrating Eq. (|39|) over (dX) l: we obtain: 



da d (Y i\ 

with the current of probability given by: 



(45) 



4ml \ X 



[(R* + (X) + R*_(X)) (R* + (X) + R*_(X))] g 
F a 2 I 1 \ dg 



2ml \X / dX 



(46) 



In the two relevant limits, R*i_(X) + R*_(X) takes the 
form: 



R* + {X) + R*_(X) = 
lQvml exp(— n X) A = 
R+(X) +R*_{X) = 
(4\ 2 m^/3K)exp(-K X X) Am > 



(47) 
(48) 



and in the limit A = we recover Eq. (4-4) of Ref. p7| . 

Important works on the solution of Eq. ( ff5[ ) , which also 
go beyond the factorization approximation by taking into 
account correlations of consecutive domains, are given in 



a series of papers by Nagai and Kawasaki 1 30 1 . Here, we 
will follow Ref. |27j and the first of the papers cited in 
Ref. 

The time dependence of the density of kinks n(t) —or 
alternatively of the average kink-kink distance: X(t) = 
(X) = l/n(t)— is studied by assuming that at large times 
X represents the only relevant scale in the problem, and 
therefore g(X;t) satisfies the scaling expression: 



g(X;t) = n(t)g(X/X) 



(49) 



For example, for a Dirac-delta distribution (all the do- 
mains have the same size) g(s) — 5(s— 1), and for a Pois- 
son distribution (randomly distributed kinks) g(s) = e~ s . 

Secondly, we will use a steady-state approximation |2jj ] 
according to which the distribution g(X; t) does not de- 
pend on time, on scales sufficiently small with respect to 
X(t): more precisely, on scales X < X* . This means 
that the motion of a couple of kinks at distance smaller 
than X* is essentially independent on the position of all 
the other kinks. Because of the scaling hypothesis, it 
must result: X* = X/a, with a constant. 

The temporal variation of n(t) is determined by the 
number of kink-kink annihilations per unit time and unit 



length. Since each annihilation makes two kinks disap- 
pear, we have: 



hit) = 2n{t)j{X = Q;i) = 2n(t)j(X*;t) 



(50) 



where the second relation derives from the fact that dtg = 
implies dxj = 0. 

By approximating (f(X)} with f(X), (/ is a generic 
function), and by neglecting R^(X) with respect to 
R%(X*), we finally obtain the following expression for 
the current in X : 

j(X*;t)= -^L[R* + (X*)+R*_(X*)]g(X*) 



F a 2 1 dg_ 
2ml X dX 



1. Deterministic regime 



If the noise term is negligible 
a 1 



3(X*;t) 



4ml X 



[RUX*) + R*_(X*)]g(X*) (51) 



Let us consider separately the two limiting cases. When 
A = 0, by using Eq. (|47j), the current writes: 

j(X*;t) = -4avg(l/a)n 2 (t)exp(-n /an) (52) 



and Eq. (|5C|) becomes: 

h(t) — —8avg(l/a)n 3 (t)ex-p(—Ko/an) (53) 
whose solution gives, at large times: 



X(t) ~ (a/K ) ln(t/ti 



8g(l/a) 



5> (A - 0) 



(54) 



In the opposite limit of a strong symmetry breaking 
(Amp V vK), a similar calculation gives: 



X(t) ~ (a/K X )ln(t/t 2 



t 2 = 



3ea 2 



4g(l/a) 



K 



(Am > 



(55) 



We therefore obtain that t\ ~ t 2 ~ t*, where t* was 
defined in Sec. Ill as the time necessary for the developing 
of the linear instability of the flat surface. So, the time 
scale for the logarithmic coarsening doesn't depend on A, 
but the length scale does, since it depends on the width 
of the (largest) domain wall. 

We can ask what is the meaning of the a-dependence in 
Eqs. ( p4}]55| ). As pointed out by Nagai and Kawasaki |30| , 
since a In t = \nt a the parameter a should have some 
"universal" value. In a mean-field calculation these au- 
thors find a = 1, while in a numerical solution of the 
kink equations they obtain a ~ 3.5. More rigorous cal- 
culations |3(J give a = 2.27 if domains are completely 
uncorrelated, and a = 3.56 if correlation effects between 
neighbouring domains are taken into account. 



8 



2. Noise-dominated regime 



Now the current is: 



<y* *\ F ° a2 1 d 9 
j{X ;t) = -—2-== — 

2TOq X oX 



The equation for n(t) writes: 
F a 2 



n(t) = - 
and the solution is: 



$'(V") 



(56) 



(57) 



X{t) 
X = 



X (t/t ) 1/3 
3a~g'(l/ar 1/3 



1/F a 



(58) 



So, we will have logarithmic coarsening at "small" 
times and a power-like one at later times. The cross- 
over time is determined by the relation (a/n) \n(t c /t*) — 
Xoitc/to) 1 / 3 ■ By neglecting the logarithmic dependence 
(also because t* 3> io), it is found that: 



a 

to | 

■ V 



(59) 




So, the ratio between the cross-over time in presence of 
a strong asymmetry and the cross-over time in absence 
of the A-term is approximately given by: 



(60) 



It is important to stress the cubic exponent in the pre- 
vious expression: so, even a not large value of {k-o/k\) 
gives rise to a logarithmic coarsening which proceeds for 
a much longer time, because kink interaction is stronger 
and therefore a larger t c is necessary so that noise get the 
better of the deterministic regime. 

We want to emphasize that in the noise-dominated 
regime, the actual value of a is much less relevant than in 
the deterministic regime, because of the power-low char- 
acter of the coarsening. 



VIII. DISCUSSION 

The main result of the present paper is that "coarsen- 
ing laws" don't change if the symmetry-breaking current 
jsB is put in the problem (at least, as far as a continuum 
local description is valid: see below). This is mainly due 
to the fact that the functional form of the kinks does not 
change, as shown by the exact solution we have given in 
Sec. |[y] for their profile. 

So, a first question is how general is this result if we 
modify the surface current, and therefore Eq. (J17|). A 
first obvious modification would be to replace d x A(m 2 ) = 



\mm! with a more complicated expression of the slope 
to. This corresponds to have a A depending on m; in 
fact, A = A (to 2 ) = 2A'(m 2 ). Since in the late stages 
of growth the slope is almost everywhere equal to ±?Bo, 
A is almost everywhere a constant equal to A(?tiq). Is it 
possible to simply replace A by A (to 2 ,) in the final results? 
This should not be a bad approximation, as suggested 
by the analysis of Eq. ( |l7j ) when A depends on m. In 
fact, the asymptotic behaviour of M{x) [the relevant one 
for kink interaction] and the values of k± can be found 
by linearizing the differential equation with respect to 
(too — M(x)) for a positive kink and to (—mo + M(x)) 
for a negative kink [in both cases, in the limit x — > oo]. 
Because of the linearization, only the value A (to 2 ,) enters 
in the problem and therefore determines the profile. 

In a similar way, we can take into account a possi- 
ble TO-dependence of the quantity K. In this case, such 
dependence might arise from a slope-dependent mobil- 
ity r [pl| —iiK has an equilibrium origin— or from the 
dependence on the terrace length I of the probability to 
nucleate a new terrace if K derives from nucleation 
noise. 

Let us now discuss the choice of the slope-dependent 
current: jes = vm(\ — m 2 /to 2 ,). The only features we 
require to have a phase separation process are: j'Es( m = 
0) > (to make the flat surface unstable) and jEs{mo) — 
for some finite value too [indeed, Too must be the first 
zero of Jes]- These features define the so-called model I. 

Modifications of ]es inside this model do not change 
the given picture, as suggested by the analysis of the 
stationary profile of the kink (for the sake of semplicity 
we put A = 0). If we linearize the equation: 

j E s(m)+Km"(x) =0 (61) 

with respect to e(x) = mo — m(x), we obtain: 

j' ES (m o )e(x)+Ke"(x) = (62) 

whose solution is again an exponential function. So, for 
x — > oo: m(x) = mo — eoe~ KX , with k = y [ /—j' ES (mo)/K. 
In our expression of jes (Eq. 10): j'Es( m o) = —2^ and 



k reduces to kq = \f2v/ K. 

Conversely, in model II there is no finite zero in ]es- 
This implies that the slope increases with no upper limit: 
for A = 0, as shown by Hunt et al. [^3|, the maximal slope 
M in the profile is asymptotically proportional to the 
size of the mounds: Mo(t) ~ X(t). Since the potential 
energy U(m) [[/'(to) = — Jbs(to)] has no minima, it is 
no more possible to define domains and domain-walls, 
i.e. kinks. 

Concerning the time dependence of coarsening, the 
only existing numerical results are the ones found by 
Hunt et al. p3| . According to their simulations (in 
presence of noise), X{t) « t n with n ~ 0.22: a fairly 
small value |52|. No (rigorous) theoretical derivation of 
n is available at the moment. Some scaling arguments 
—applicable to noiseless growth— can be found in Rost 







and Krug (3^] and in Golubovic J34]: the former give 
n < 1/4 while the latter gives the equality n = 1/4 j35| . 

A final question we want to face now is how narrow 
kinks A actually are. In the limit Atoo 3> \pvK, from 
(21) we have: n + = Xtjiq/K and /t_ = 2^/Atoo. A 
e inspection shows that [A] = [K] 

A 



Eq. 
simp 



= length" 
K^F 



time . 
%. This 



Previous evaluations suggest ]6jl6 
expression for is surely wrong, if thermal detachment 
plays an important role. Conversely, if A and K —or, 
more precisely, A(?t1q) and K(rn^)— are of the same or- 
der of magnitude, we obtain re + w too- This means that 
the width of the positive kink (— l/«+) is nothing but 
the inverse of the value of the constant slope in the sur- 
face profile: so, if too is determined by the symmetry of 
the crystal lattice, too — 1/a and the positive kink is as 
narrow as a lattice constant! In this case, our descrip- 
tion would break down, because the regions of positive 
curvature in the surface profile would correspond to a 
discontinuity of the slope, i.e. to angular points, which 
are not compatible with a local continuum equation p6[ . 



Muffins term (Km"(x)) is replaced by a higher order one 
(Km"" (x)) and the resulting equation is studied numer- 
ically. Clearly, in two dimensions analytical treatments 
are much more difficult; anyway, neither a numerical so- 
lution of the model studied in the present paper is avail- 
able at the moment. One reason is that in two dimen- 
sions, even the model without jsb is not yet fully under- 
stood, since the evolution equation for m(x, t) ||37|| is no 
more equivalent to the Cahn-Hilliard equation. 
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APPENDIX A: LANGEVIN EQUATIONS FOR 
THE KINKS 



IX. CONCLUSIONS 



1. Absence of noise 



The kink picture not only has allowed to find the coars- 
ening law in presence of the symmetry breaking term, but 
it has also given a qualitative description of the dynam- 
ics which allows a better comprehension of the evolution 
of the system: the widening or the narrowing of a kink; 
the consequent different velocities of kinks A and B; the 
conservation of the order parameter seen as a kinemati- 
cal constraint on kink movement; the difference between 
the "real" force acting on a kink and the "effective" force 
felt by the kink, because of such constraint. 

In this respect, the most important consequence of the 
breaking of symmetry is that negative kinks feel an ef- 
fective repulsive interaction with the nn kinks (but at- 
tractive with the nnn ones). It is important to stress 
this point because coarsening is the result of a global at- 
traction between kinks: if kinks repelled each other, the 
configuration with the Xi all equal would be stable. 

Finally, the kink picture has provided the condition of 
applicability of the local theory: 

— = b \ » a 63 

TOoA(TOg) 

If this relation is not fulfilled, a different method to study 
coarsening should be used. In Ref. [|| we showed that 
in this case the evolution of the surface is governed by 
a nonlocal current; alternatively, we can keep a local de- 
scription, but we must add a singular term to the current 
j, and couple the Langevin equation: dtz(x,t) — —ad x j 
with specific evolution equations for the angular points. 
It would be clearly interesting to check if a different coars- 
ening process may arise from an "angular point" picture. 

To our knowledge, the current (|0|) has not been for- 
merly studied. The closest model is the one consid- 
ered by Stroscio et al. Jl8| in two dimensions, where the 



The starting point is the following multi-kink expan- 
sion: 

m(x,t) = Mi(x,t) + ^2[Mj(x,t) - Mj(-oo)] 

j>i 

+ Y,[M j (x,t)-M j (oo)] (Al) 

(A2) 



Mi(x,t) + Snii 



which gives rise, once replaced in (|22|), to: 

a_1 E 



dvj 



= D x j[m(x,t)] (A3) 



Mj depends on x and t through the combination (x — 
Vjt) and Mj is the derivation with respect to all this 
argument. We will also use the notation: d x M to mean 
the same kind of derivation. The single kink profile is 
found by simply dropping the sum Y) j and the term in 



v in Eq. (A3) 



- a-^M'j = DvjlMj] (A4) 

It will be useful to consider, together with Mj, also 
its spatial derivative Mj which is localized around x = 
Xj. We define also Mj through the relation: Mj(ai) = 
D x Mj(—x). They satisfy the relations: 

-a^v^M-ix) = 

T> x [Kd x - U"{M 3 ) + AMj + AM i d x ]Mj(a;) (A5) 
a^v^M'jix) = 

[Kd 2 x ~ U"(Mj) + AM' + XMjdxpxM'Ax) (A6) 
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Now, let us multiply (|A3|) by M-(x) and integrate on 



x. By defining Svj — i 



in, we can write: 



dxM<D x j[m] + a 



dxM-M'j 



(A7) 



The next step is to replace m(x, t) = Mi(x) + 5m, in 
the current j [to] . The definition of the nonlinear part of 
the potential [/(to) (or equivalently of the current U'(m)) 
is self-explanatory. 

j[m] = Km" — [/'(to) + Atoto' 

= KM" + KSm'l - U'(Mi + Sm,) 

+\{M l + 5m i ){M' l +8m' i ) 

= KM'l + KSm'l - U'[Mi) - U" '(M^Sm, - U' NL<i 

■ + XSrriiMl + XMiSm^ + XSmrfm^ (A8) 

The th ree terms which do not depend on Smi, once 
used Eq. (A4) cancel the term j = i in the last summation 
of Eq. QA7D; 8m" is simply written as Ylj&Mj, whilst 
all the other terms for the mom ent keep unchanged. So, 
the right hand side of Eq. (A7) rewrites: 

RHS\(A7) = 

J dxM[{x) \ YjKV x M'l + o-^Mj] 
-V x [U' NL ^ + U"{M i )8m i ] 

+XD x [Sm l M- + MiSm'i + Sm^m-]} (A9) 
Let us consider separately some terms: 

dxM-(x) ^2 K~D X M" = J dxKSmid 2 x D x M{ 



- / dxM-(x)D x [U"(Mi)Smi] = - J dxU"(Mi)S mi D x M^ 

X J dxM-D^SmiM'i + Sm'iMi] = 

X J dx[8m i Ml + 2M i 8m' i }T) x M' i 

+A J dxSmi[M- + Mid x ]D x Ml 
Eq. (AG) therefore rewrites: 

RHS\(A7) = 

dx8m l [K& 2 x - U"{Mi) + XM- + XMid x ]D x M- 



f dxMlJ2 a ~ lv ° M l 



This way, Eq. (|A7|) takes the form (LHS ^left-hand- 
side): LHS\U) = (A) + (B) + (C). 



By using Eq. ( 

(A) = a 



1 J dxMCZy^M^ (All) 



which can be summed to (£?), giving: 

(A) + (B) = a- 1 / dxMlJ2(v°j ~ «?)Mj 



(A12) 



and subtracting LH S\ : 

(A) + (B)-LHS\^ 
that is to say: 



_1 E 



dxMl[D x U'js 



NL. 



-XD x (SrmMl + 2M,8m' + Smrfm'J} 
(Ci) + (C 2 ) + (C 3 ) + (Gi) 



(A13) 



In the previous equation we have used the following scalar 
product: 



(R, S) = 



+oc 



dxR(x)S(x) 



(A14) 



The three terms in square brackets on the right-hand-side 
[(C 2 ) + (C 3 ) + (C4)] represent the effect of the symmetry 
breaking current. 

By integrating by parts and by using the definition of 
M': 



(d) = (Ml(-x),U' NL A 



(A15) 



+ / dxM{ [-D x U' NL i + XD x (8m,iMl + 2Mi8m' + Jto^to-)] 
= {A) + (B) + (C) (A10) 



If we define the function G(x,y) = U'(x + y) — U'(x) — 
yU"(x), then U' NL i — G(Mi,8rrii). In the following, we 

will also make use of the function: U(x, y) = U(x + y) — 
U(x) — yU'(x). It is obvious that G(x,y) = d x U(x,y). 

We observe that: 1) G(Mj,0) = 0; 11) G may be writ- 
ten as a Taylor expansion whose generic term contains 
{8mi) n ; 111) G is not linear in Srrii, but if we use the bi- 
nary interaction approximation, it is indeed linear. This 
approximation corresponds to write: 

(Sm t r « Y^M, - M^-oo)]" + YjMt - M 3 .(oo)]» 

j>i j<i 

In this approximation, we obtain: 

(Ci) = YjJ dxM i( x i - x)G(Mi(x - Xi),Mj - Mj(-oo)) 

j>i 

Y J dxM-(xi - x)G{M t {x - Xi), Mj - M,-(oo)) 
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We must observe that M[ is not vanishing only when 
x « Xi\ furthermore, [Mj — Mj(±oo)] goes to zero when 
(+) x > Xj or (— ) x < Xj. On the basis of these consid- 
erations, it is possible to write: 

r+ao 

(Ci) = J2 / dxM[( Xl - x)G(M l (x - Xi),AMj) 
j>i ■'• 

+ 22 dxM-(x l - x)G(M l (x - Xi), -AM,-) 

• ^ • J — OO 



where AMj ee Mj(oo) - Mj(-oo). 
Since M-(x) is an even function of x: 



because the second one would be exactly zero). We give 
here only the results. 



dxM[(xi — x)G(Mi(x — Xi), const) = 
U(Mi{x — Xj), const) 

and (Ci) can be written as: 
(GO = 5^[&(Mi(oo), AM;) - U(Mi(xj - Xi),AMj)] 

j>i 

+ Y}U(M t ( X] - Xi ),-AMj) - U(Mi{-oo),-AMj)} 



3<i 



At the first order in the small quantities [Mi(xj — Xi) — 
Mj(±oo)] (± resp. for j > i and j < i), we have: 



(Ci) = - Yl\. M i{xj ~ Xi) - M ( (oo)]G(M i (oo), AMj ) 

j>i 

+ Y^[Mi(xj - x t ) - M l (-oo)]G(M l (-oo), — AMj) (A16) 

In the following, we will restrict ourselves to nearest- 
neighbour kinks interactio n, an d therefore only the terms 
j = i ± 1 will survive in (Alf ). If we also use the fact 
that: 

G(Mi (±00), ±AM i± i) = TAMi±iC/"(Mi(±oo)) (A17) 

we obtain the following final expression: 

(d) = [Afi(x i+ i - x t ) - Mi{oo)]AM i+1 U"{Mi(oo)) 
+[M i (x i _ 1 - xi) - Mii-oo^AMi^WiMii-oo)) (A18) 

The procedure to follow for the treatment of the other 
terms (Ci) is similar. In poor words, if R(x) and S(x) 
are functions which are localized, resp. in x\ and X2, we 
make the approximation: 

R(x — xi)S(x — X2) ~ R(x — xi)S(xi — X2) 
+ R(x 2 — xi)S(x — X2) 

and then we retain only the term corresponding to the 
function decreasing more rapidly (for example, if R(x) 
was a Dirac-delta, only the first term would be retained, 



(G 2 ) = \{[M i+1 (x i+ i - Xi) - M i+ i(oo)] 



+ OO 



- [Mi-x(xi - Xi-x) - Mi_i(oo)]} / dx(M^) 2 

J — OC 

(C s ) = 

(GO = -^(AAf J+ i) 2 M;(x 4+1 - Xi ) 
+ ^(AMj_i) 2 M/(sj — Xi_i) 

The expression (G3) = means that such term is always 
of higher order than the others. 

From now on, notation must take into account explic- 
itly the existence of two different classes of kinks. By 
using the following results: 

Mp(x) — Mp(oo) ~ — /32mo exp(— Kpx) when x — ► 00 

Rp(x) = exp(— K/32) 

AM 4 = /32m 

t/"(±m ) = 22/ 

+00 2 

dx[M^(x)] 2 = -m^Kji 

-do " 

it is straightforward to write: 

(GO = %v m l\Rp(X,{) - RpiXi.t)} 

(G 2 ) = P(4/3)m 3 QKf3 X[R- (X l ) - R-piXi-j)] (A19) 

(G 4 ) = -/34™^ /3 A[i? /3 (A l ) - ^(Xi-i)] 



We can now put together Eq. ( |A19| ) with the previous 
ones. A further approximation is to neglect the "defor- 
mation" of the kink-profile, due to its velocity, and to 
suppose that kinks are immobile in absence of interac- 
tions. This way, we obtain: 

a' 1 E v i(Ml, M'j) = (G x ) + (G 2 ) + (G 4 ) (A20) 
3 

where the LHS can be further developed: 



a-^Vjj dxMl(x)M^(x) = 
3 

a^Y^Vj J dxB x M^(x)D- 1 M^x) = 
3 

a^Y^vj J dxMl(x)T>- 1 M' j (x) (A21) 



We have therefore to determine the inverse of the op- 
erator D x . By following Kawasaki and Ohta pa: 
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D- 1 A(x,t)=~ J dx'\x- x'\A{x',t) (A22) 

The "integration constants" appearing when the operator 
is inverted, are shown to be irrelevant for the kink 
dynamics (Ref. p7| 



By applying Eq. (A22) to Eq. (A21) 



dxM'^D-'M'ix) 



--\xi-xj\AMiAMj 



dxdx'M<(x)\x-x'\M' j (x') (A23) 

-|.r, -.v, I I <i.r,lx'Mi(x-Xi)Mj(x-Xj) (A24) 

(A25) 

So, Eq. ( |A13| ) finally writes: 

3 

(A26) 

2. The effect of noise 

The term of noise SF(x, t) in Eq. (||) corresponds to 
a term rj(x,t) — d x SF(x,t) on the right-hand-side of 
Eq. ([A3|). To see how it affects the kink movement, it 
must be mul tiplie d by M^ix) and integrated on x. Since 
the LH S of ( |A2fj| ) indeed corresponds to minus the LHS 
of (A3), if we call r)j(t ) the noise term to be added to 
(Ci) + (C 2 ) + (C 4 ) in (|A26|) , it will result: 



r]i{t) = - J dxMl(x)ri(x,t) = J dxd x M-(x) ■ 5F(x,t) 

(A27) 

The following properties are found [^7j : 
(%(*)) = and 



(0) =2F 5(i-«') / dxM'(a;)D x Mj(a;) 

= -4m§F (-l) i_ -'>i - a:jl5(* - t') (A28) 

To derive the spatial correlation between noise, we have 
used the definition of M- and inverted the operator 
D x . Finally, we have used the fact that AMiAMj = 
4mg(— , a rel ation which can be used also for the 
LHS of Eq. ( A26| ). So, we obtain the following system 
of coupled Langevin equations: 



-2a 1 m%y](—iy 3 \x i -x j \xj = 

3 

(C 1 ) + (C 2 ) + (C 4 ) + r?i(t) 



(A29) 
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